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ABSTRACT 

We perform a detailed 2-dimensional weak gravitational lensing analysis of the nearby (z = 0.058) galaxy 
cluster Abell 3128 using deep ugrz imaging from the Dark Energy Camera (DECam). We have designed a 
pipeline to remove instrumental artifacts from DECam images and stack multiple dithered observations without 
inducing a spurious ellipticity signal. We develop a new technique to characterize the spatial variation of 
the PSE which enables us to circularize the field to better than 0.5% and thereby extract the intrinsic galaxy 
ellipticities. By fitting photometric redshifts to sources in the observation, we are able to select a sample of 
background galaxies for weak lensing analysis free from low-redshift contaminants. Photometric redshifts 
are also used to select a high-redshift galaxy subsample, with which we successfully isolate the signal from 
an interloping z = 0.44 cluster. We estimate the total mass of Abell 3128 by fitting the tangential ellipticity 
of background galaxies with the weak lensing shear profile of an NEW halo, and also perform NEW fits to 
substructures detected in the 2-D mass maps of the cluster. This study yields one of the highest resolution mass 
maps of a low-z cluster to date, and is the first step in a larger effort to characterize the redshift evolution of 
mass substructures in clusters. 

Subject headings: galaxies: clusters: general - galaxies: clusters: individual (Abell 3128) - gravitational 
lensing: weak - techniques: image processing 


1. INTRODUCTION 

Cosmological perturbation theory provides a framework in 
which cold dark matter organizes itself hierarchically, first 
collapsing into small structures which can overcome cosmo¬ 
logical expansion and then continuing to merge into increas¬ 
ingly large halos. Because small collapsed objects often sur¬ 
vive accretion onto a larger system to become sub-halos of 
their host, the hierarchical structure formation paradigm pre¬ 
dicts that dark matter halos should be ric h in mass substruc¬ 
tures ( [Klypin et al.|2011[|Gao et al.|2Q12 ). 

The amount of mass substructure that we observe should 
increase with redshift, particularly for cluster-sized halos: 
mergers of galaxy- and group-size halos are more common 
in the early Universe, and the resulting clust ers have long dy¬ 


namical relaxation times ( |Gao et ap|2004| ). At higher red 


shifts, then, we expect to observe an increasingly high fraction 
of the total cluster mass locked up in sub-regions of localized 
mass enhancement. An observational study of cluster sub¬ 
structure and its cosmic evolution would probe the assembly 
history of cluster-sized halos and test the CDM paradigm on 
sub-megaparsec scales. Characterizing substructure in clus¬ 
ters also has important implications for understanding the role 
of the mass environment in the evolution of member galaxies. 
Correlating sub-halo locations and, e.g., star formation rates 
would reveal the effect of local mass environment (distinct 
from the larger-scale mass distribution) on galaxy properties. 

Studies of cluster substructure are already underway. Most 
notably, X-ray data have been used to obtain cluster mass 
functions through the proxies of cluster gas emissivity and 
temperature. However, these proxies are related to mass by 
scaling relations that rely on assumptions like the hydrostatic 
equilibrium of the intra-cluster medium. Accretion-induced 
heating of cluster gas, as well as merger-induced bulk and tur¬ 
bulent motions, violate the assumption of hydrostatic equilib¬ 
rium at the substructure level, at which the cluster is dynami- 
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cally active. Moreover, the error induced by the assumption is 
greater in the outskirts of higher redshift clusters, where merg- 
ers are more frequent and clusters ar e accreting more rapidly 
( |Lau et al.|201 3 [[Nielson et al.|2014| ). Hence, an accurate ob¬ 
servational study of cluster substructure requires an analysis 
technique insensitive to the dynamical state of a cluster. 

Because of its freedom from assumptions about baryonic 
physics, weak gravitational lensing (WL) has become a stan¬ 
dard tool for measuring mass concentrations in the Universe. 
Multiple observations show that individual dark matter sub¬ 
structures within a cluster are c apable of producing their own 


detectable weak lensing shear (|Qkabe & Umetsu 


Se- 


__ _ Tg snear ( _ _ __ 

|hgal et aL]|2Q08 Huwe|MT^ . However, in the weak lens- 
mg regime, the distortion of background galaxies induced by 
the intervening cluster is much smaller than the intrinsic un¬ 
certainty of galaxy shape measurement. To overcome this 
challenge, obtaining the angular resolution needed to iden¬ 
tify cluster substructures, requires large numbers of resolved 
background objects. We achieve this in our observations by 
taking deep, wide-field images using the Dark Energy Cam¬ 
era (DECam) on the Blanco 4-m telescope at Cerro Tololo 
International Observatory. 

In this paper, we describe the analysis pipeline we devel¬ 
oped to make weak lensing measurements on DECam data 
and present the results of our first substructure study on a 
nearby (z = 0.06) cluster, Abell 3128, which is one of the 
hig hest resolution mass maps of a z < 0.1 cluster to date (see 
also lQkabe et al.|2014| ). We chose to begin our investigation 
of cluster substructure with Abell 3128 for several reasons. 
Eirst, weak lensing analysis requires clusters to be massive 
enough that their (obviously less massive) substructures can 
be detected, and A3128 is one of the most massive clusters in 
the Local Volume without a published convergence map. Sec¬ 
ond, the large projected size of low-z clusters (and their sub¬ 
structure) means that their distortions will be coherent over a 
large swath of the sky. This greatly increases the number of 
background galaxies from which to measure the WL signal, 
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and makes substructure easier to detect at low redshifts. Fi¬ 
nally, Abell 3128 has a complex morphology that has been 


well studied at radio, optical and X-ray wavelengths ( [Rose et 
al.|2002t Werner et al.|2007| ), enabling the comparison of our 
WL substructure analysis to other techniques. 

Our observations of Abell 3128 and reduction method are 
described in §3 and §4, respectively. In §4 we also dis¬ 
cuss our weak lensing analysis and substructure identification 
schemes, whose results we present in §5. In §6 we explore 
the significance of our findings. Finally, we summarize our 
conclusions in §7. 

2. WEAK GRAVITATIONAL LENSING 

A thin gravitational lens deflects and distortsthe images of 
background sources. It is customary to describe the mass dis¬ 
tribution of the lens in terms of its surface mass density S or 
its dimensionless convergence. 


K = 


( 1 ) 


where the critical surface mass density l^crit of the lens is de¬ 
fined as 

Ds 

=-(2) 

AnGDiDis 

The quantities Ds, Di and Dis are the angular diameter dis¬ 
tances to the background source (galaxy), the lensing cluster, 
and between the source (galaxy) and lensing cluster, respec¬ 
tively. 

Lenses with a: > 1 produce the dramatic arcs, Eins tein rings, 
and m ultiple images of strong gravitational lensing ( |Coe et ah 


|2Q10| ). Areas of a: < 1 define the weak lensing (WL) regime, 
in which the distortion of background images produced by the 
lensing cluster is much smaller than the background images 
themselves. Because a: > 1 typically only in the dense core 
of a cluster, to study substructure over a cluster’s entire virial 
region requires WL analyses. 

The convergence k produces the isotropic magnification of a 
background source, and as such cannot be measured directly 
from an image without some prior knowledge of the source 
size. In the case of weak lensing, however, convergence can 
be recovered from the shear y of background source images 
caused by the tidal forces of the lens’ gravitational field. In 
particular, convergence is related to shear through the applica¬ 
tion of the inverse 2-D Laplacian to the 2-D gravitatio nal po¬ 
tential of the lens , an i ntegral transform fir st derived in |Kaiser| 
& Squiresl (119931) andlKaiser et al. I (119951), with variants pub- 

|l4hlnian ^ ' * 


Fahlman et al.|(|1994|) and 

In the WL regime {k << \,y « 


Fischer & Tyson|dI997| ). 

, the background galaxy 


images experience a curl-free stretching in the direction tan¬ 
gential to the line-of-sight from the lens. Now, galaxies are 
elliptical objects with a measured ellipticity e = 1 - (bio) 
tilted at some position angle 6 with respect to the image axes. 
Galaxy shapes are frequently decomposed into ellipticity mo¬ 
ments: ei = ^cos(2^) is the projection of a galaxy’s ellipse 
onto the image x and y axes, and e 2 - e sin(2^) is its pro¬ 
jection onto the lines y = x and y = -x. With this in mind, 
the the shear induced by a lens can be written in terms of the 
galaxies’ shapes as 


r ^ ^tan = -(e\ cos(20) -r ^2 sin(20)). 


(3) 


In this so-called tangential ellipticity, (p is the angle from a 
fiducial lens center to the galaxy measured counterclockwise 
from north. In other words, the factors el and e2 characterize 


a galaxy’s shape and position angle relative to the image axes, 
and the factors of sin(20) and cos(20) rotate the galaxy’s e\ 
and e 2 into the line tangent to the radial extending from the 
chosen lens center. In the absence of a gravitational lens, 
the average tangential ellipticity <^tan) should vanish when 
considered over many background galaxies with no intrinsic 
alignment. Hence, the <^tan) measured at a point in the image 
is an unbiased estimator of the WL shear. By measuring the a 
systematic deviation from zero average ellipticity with a sam¬ 
ple of galaxies widely separated from each other in redshift 
space, we may reconstruct a cluster’s convergence. If in ad¬ 
dition the redshifts of the cluster and background galaxies are 
known, the cluster’s surface mass density ma y be recovered. 
For a comprehensive treatment, see reviews by|Bartelmann & 
|Schneide^ ( |2001| ) and | Wittnian] ( |2002a| ) . 


In this study, we identify and characterize mass peaks us¬ 


ing the ape rture mass statistic Map first introduced by Schnei 
der| ( [I996] ). Measured some angular distance 6 away from the 


cluster center. Map is given by the convolution of the conver¬ 
gence K with an aperture mass filter f/(|^o “ ^D- 


MapiOo) = I J (feK{0)U{\9o - 01). (4) 

The aperture mass filter f/(|^o - ^1) smooths the convergence 
over some characteristic aperture Oq. By design, the Map is 
a local measurement involving only the shear from galaxies 
within an angle of the center at position^; the filter is “com¬ 
pensated” so that its first order moment vanishes on scales 
larger than the aperture size. 

If the aperture mass filter t/(|^o - ^1) in Equationj^is trans¬ 
formed as 

2(100 - 0|) = I 0V(|0o - - u(0), (5) 

we can replace k in the aperture mass statistic with the tan¬ 
gential ellipticity of background galaxies ^tan- For a discrete 
dataset of background sources, the aperture mass thus has the 
form 

MapiOo) = ^ X ‘'r(0)2(|0o - 01), (6) 


where the sum is taken over all galaxies in the observation 
( |Schneider|| 1 996 ) . We apply Equation to our observations 
to build WL convergence maps. 


. convergence maps. 

A variety of aperture mass filters exist which can be used 
in Equatioij^ exist, but the one best suited to our sea rch for 
substructure was introduced by|Schirmer et al.|(|2004|) as part 
of the GaBoDS survery. The Schirmer filter was originally 
designed to pick out shear signal from clusters embedded in 
large-scale structure, and is given by 


Q(x) = 


1 


tanh(x/xc) 


(1 -r e^-^^ + tiRUxIxc) ’ 


(7) 


where Rs is the Schirmer filter size and x = r/Rs is a scaled 
distance between the cluster cent er and the point in consid¬ 
eration ( [Hetterscheidt et al!]|2005| ). To optimize the filter for 
detection of NEW shear profiles, the parameters in Equation 
[7]are tuned to a = 6,b = 150,c = 47, J = 50 and x^ = 0.1. 
By tuning the size of the Schirmer filter (Rs) in the aper¬ 
ture mass statistic, we can discern both the main cluster sig¬ 
nal and its substructures while simultaneously characterizing 
their respective scales: noting that the Schirmer filter weights 
peak sharply at a value of XcRs , the structures identified have 














































Mass Substructure in Abell 3128 


3 


size ~ O.IO/?^. Since the Schirmer filter is not monotoni- 
cally decreasing, it is difficult to assign it with a Gaussian- 
type FWHM. Instead, a smoothing length can be obtained by 
computing the radius which encompasses 50% of the filter’s 
weight. For the form of the Schirmer filter given above, this 
radius is 0.121/?^, equivalent to a Gaussian with standard de¬ 
viation cr = 0.18 or a FWHM of 0.3. 

We note finally that, in analogy with electromagnetism, the 
curl-free tangential ellipticity defined in Equation is some¬ 
times called the E-mode WL signal. A complementary, curl¬ 
like ellipticity referred to as the B-mode signal is obtained by 
rotating Equation [^through nj A radians: 

ec = €2 cos(20) - ei sin(20). (8) 

Since gravitational lensing creates no B-mode signal, the re¬ 
placement of ^tan with ec is frequently used as a statistical 
control. 


3. OBSERVATIONS 


The DECam imager consists of 62 2048 x 4096-pixel sci 
ence CCDs (520 megapixels total) arranged in a hexagonjand 
captures 3 square degre es (2.2 square degrees wide field) at 


0.265”/pixel resolution ( |DePoy et al.||2008] [Elaugher et al. 
|2012| ). The camera’s wide field of view allows us to image 
the entire virial region of even a low-redshift cluster in a sin¬ 
gle pointing, making it efficient for our study of cluster sub¬ 
structure. 

Observations of A3128 were made over eleven days from 
8th-24th November 2012 in the ugrzY filter set by Dara Nor¬ 
man and the DECam science verification team as part of that 
instrument’s science verification program. To ensure sky cov¬ 
erage in the gaps between science CCDs, the telescope was 
dithered in a “center -r rectangle” pattern. The dithers are 
large enough to overlap adjacent CCDs by several hundred 
pixels, providing more uniform depth at the chip edges and al¬ 
lowing construction of a catalog covering the entire 1.5° x 1.5° 
uniformly. The exposure time of each pointing varied by fil¬ 
ter: 720 s in u, 600 seconds in g, 300 seconds in r and 240 
seconds each in Y and z. The final exposure times across the 
field were 10,800 seconds in u, 3600 seconds in g, 5400 sec¬ 
onds in r, 2630 seconds in Y and 2160 seconds in z, with at 
least two complete dithers in each band. The mean seeing in 
r was 0.94”, and after calibration of source number counts vs. 
magnitude against the Subaru-COSMOS catalog ( |TaniguchT| 
et al.|2007] ), the observations have a 50% completeness depth 


of m = 24.97 'm u, m = 24.76 in g, m = 25.62 in r and 
m = 23.85 in z. We note that the Y and z filter profiles overlap 
considerably, such that the narrower Y essentially just covers 
the longer-wavelength portion of the z filter. Imaging in these 
two filters is redundant for the purposes of our analysis, and so 
we make no use of the Y band data beyond making a stacked 
image. 

The CTIO-l-DECam system’s sensitivity is greatest in r 
band, and this filter also optimizes the balance between high 
background galaxy luminosity and reasonably low sky noise 
(both of which increase with increasing wavelength). Follow¬ 
ing the successful observi ng strategy of the Deep Lens Sur¬ 
vey (|Wittman et al.|2002bl), we observed A3128 in r when the 
seeing FWHM reached < 1.0” and in ugzY otherwise. Ac¬ 
cordingly, the r-band imaging has uniformly good resolution 
as well as a greater depth than the imaging in other bands. 


^ Note that as of the time the A3128 data was collected, one of the CCDs 
at the southern edge of the array (N30) was non-functional. 


Shear measurements are thus made exclusively in the r band, 
while other filters are use d to p rovide color information for 
photometric redshifts (see |4.2| ). 

4. ANALYSIS 
4.1. Image Processing 

Reconstructing a two-dimensional mass maps from galaxy 
shapes is an involved procedure. The intrinsic galaxy elliptic- 
ities are ~ 30 times larger than the distortions we are trying 
to measure, and a priori it is impossible to disentangle WL- 
induced shear from the shape of a single galaxy. Moreover, 
anisotropy in the PSF field shears incoming light from galax¬ 
ies and obscures our weak lensing signal. For a camera as 
large as DECam, which reaches the edge of the focal plane of 
its telescope, this effect is substantial. The number and size of 
DECam’s CCDs also makes removing instrumental artifacts 
from observations a technical challenge. In the following sec¬ 
tion, we list the image processing steps undertaken to over¬ 
come these difficulties and measure the mass substructure of 
A3128. Although a community reduction pipeline now exists, 
we developed our own independent image processing pipeline 
as part of the science verification program for DECam. 


4.1.1. CCD-Level Reduction 

The CCD-level image reduction applied to each exposure 
in the dataset includes the standard complement of overscan 
subtraction and trimming, bias subtraction, and dome fiat 
field correction. These tasks were accomplished with the 
MSCRED package in IRAI0 We apply to CCD images an 
empirically-determined correction for the crosstalk that oc¬ 
curs as neighboring amplifiers are read out in parallel. A “tree 
ring” pattern of concentric circles of light and dark pixels ap¬ 
pears in all DECam object and fiat field exposures; these are 
not an artifact of gain variations on the chips, but actually 
represent the physical shifting of charge between pixel wells. 
Tree rings in object exposures are successfully camouflaged 
by the fiat fielding step, although fiat fielding away the tree 
rings in object exposures is tantamount to turning an astro- 
metric error into a photometric one. Given the tiny amplitude 
of the tree rings (~ 0.2% of pixel flux value), the error intro¬ 
duced is dwarfed by the m > 0.03 photometric uncertainty of 
the images. For more details r egarding this and ot her CCD 
artifacts peculiar to DECam, see [Plazas et al. ( [2014 ). 

To mitigate the > 100" pointing errors in DECam sci¬ 
ence verification data, objects in the observation are matched 
against a list of reference celestial coordinates in the USNO- 
B catalog. We fit a linear relation, which may include a zero 
point shift, scale change, and axis rotation, between the ob¬ 
served positions and the reference coordinates on both coor¬ 
dinate axes. The fit is used to update the image world coordi¬ 
nate system so that it is registered to the reference coordinate 
system defined by USNO-B. 


4.1.2. PSF Correction 

Because the DECam CCD array is so large, the point spread 
function (PSF) has significant and spatially varying contribu¬ 
tions both from the curvature of the focal plane and the Blanco 
4-m optics. Such anisotropies in the point spread function can 
induce spurious shear signal, and so the accuracy of our mass 
maps relies on extremely precise characterization of the DE¬ 
Cam PSF. Distortions in the PSF field of an image can be 

^ IRAF is distributed by the National Optical Astronomy Observatory 
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traced by systematic variations in the shapes of its stars since 
these are point sources and should appear perfectly round in 
an isotropic PSF field. Consequently, the first step in circular¬ 
izing the PSF is the identification of stars in the observation. 
Next, polynomials are fit to the spatial variation of the follow¬ 
ing combinations of second-order intensity moments: I^x-Iyy, 
hx + lyy and Ixy. Finally, the intensity moment fits are used to 
derive a PSF circularization kernel, w hich is conv olved with 
image pixels in the stacking stage (cf. 4.1.3| ) as in Bernstein 
|& Jarvis] ( |2002| ). 

The default procedure is to go through these steps for ev¬ 
ery CCD of every exposure in the dataset, and for most ap¬ 
plications, this would yield a sufficiently circular PSF. How¬ 
ever, the DECam PSF field is severely under-sampled by un¬ 
saturated stars in any single exposure; the CCD chips them¬ 
selves are large, and A3128 is at high galactic latitude. Conse¬ 
quently, applying the standard PSF circularization technique 
to the DECam A3128 data only lowers the mean stellar ellip- 
ticity from 5% to 1%, which is still high enough to affect the 
WL shear signal. 

To improve the PSF modelling, we combine stellar cata¬ 
logs from sequential exposures on the same CCD chip and 
use these “super-catalogs” of stars to fit 2-D polynomials to 
the spatial variation of intensity moments. High-order poly¬ 
nomial terms (fourth- and fifth-order) of the intensity moment 
fits capture effects like focal plane curvature and tend to be 
stable over the course of contiguous dithers. Because the 
lower order terms in the intensity moment polynomials cap¬ 
ture time-dependent effects like seeing or telescope drift, they 
are usually best fit using individual CCD exposure catalogs 
rather than super-catalogs. 

The particular grouping of exposures used to build super¬ 
catalogs is also determined empirically for each CCD. Stellar 
ellipticity is minimized when the higher-order terms of the 
hx - hy polynomials are fit by super-catalogs assembled from 
a single dither of five exposures. Meanwhile, the higher-order 
terms of the Ixy and hx + lyy polynomials should be fit using 
super-catalogs assembled from as many contiguous exposures 
as possible. 

For each of the 61 functional CCDs, we determine empir¬ 
ically both the degree of the polynomial fits and which of its 
terms should be obtained through super-catalogs. We piece 
together the final forms of fits to hx - hy^ hx + hy ^xy 
from whichever combination of terms ultimately yields the 
lowest stellar ellipticities. This procedure allows for the cir¬ 
cularization of stellar PSFs to better than 0.5% (see Figure[^. 
We verified that the magnitude of stars selected does not af¬ 
fect the polynomial fitting by examining the residuals of the 
PSF fits to stars. These show no discernible trend in the range 
of magnitudes considered (16 < < 20 ). 

The effectiveness of this PSF circularization scheme may be 
quantified by constructing two-point shear correlation func¬ 
tions, defined as 


Qj = {ei(r)xej(r^6)}, (9) 

where Cf is the ith ellipticity component of an object at posi¬ 
tion r, and brackets denote an average over all pairs within a 
separation 0. A third correlation function, 

C 3 = (^i(r) X e 2 (r + 6 >) -r e 2 (r) x ^i(r -r 6)), (10) 

should be zero and is frequentl y used to test for syst ematic er¬ 
ror in PSF correction schemes ( [Massey et al.|2005] ). Star-star 
auto-correlation functions and star-galaxy cross-correlation 


functions are shown in FigurelS The star-star auto-correlation 
functions Cn and C 22 (top left panel) show a small signal on 
small scales that we attribute to some over-correction of the 
DECam PSF in the areas near stars. The star-galaxy corre¬ 
lation functions (bottom panels) show the same correlation 
on small scales, with a magnitude several times smaller than 
the galaxy-galaxy auto-correlation signal (which traces WL 
shear). The negligible value of the C 3 in both star-star and 
star-galaxy pairs confirms that the PSF circularization scheme 
introduces no major systematic error to shape measurement. 

It should be noted that even after circularization, the mea¬ 
sured ellipticity moments do not yet represent the true shapes 
of the galaxies. Both atmospheric seeing and the circulariza¬ 
tion of the image PSF make galaxies appear more round than 
they really are, and effectively dilutes the WL shear signal. 
We corr ect for this “smearing” at the catalog level as detailed 
in^ 


4.1.3. Stacking 

Once we have corrected for CCD artifacts and, in the case 
of the r-band exposures, obtained polynomial fits to the PSF, 
we proceed to stack CCD exposures into a single image. Our 
procedure for stacking closely follows the on e used in the the 
Deep Lens Survey; see jWittman et al.|2006| for full technical 
details. The steps undertaken to produce our stacked images 
are summarized here. 

1. Source Detection (yr Characterization. - Fo r each 
CCD image in th e dataset. Source Extractor ( jBertinj 
& Arnouts||1996|) is used to make a catalog of high 
SIN objects. Source Extractor also generates the sky 
background-subtracted images that will be the final in¬ 
puts to the final stack image. Subtracting the sky 
background at this stage eliminates the need to match 
sky levels at the st acking stage. At this st age, the 
ELLIPTO program ( [Bernstein & Jarvis]|2002| ) is used 
to measure the so-called “adaptive” second-order mo¬ 
ments of objects in the CCD image. Adaptive mo¬ 
ments are centrally weighted by an elliptical Gaussian, 
and their measurement is equivalent to finding the best- 
fit elliptical Gaussian for each object. Unlike SEx- 
tractor’s intensity-weighted moments, which are com¬ 
puted within some limiting isophote, adaptive moments 
do not depend on magnitude. This property makes 
adaptive moments more advantageous for g alaxy shape 
measurement and the identification of stars ([Wittman et[ 
al.|2006l >. 


2. Star Identification. - As discussed in |4. 1.2 


we use 

stars to trace out the spatial variation of the PSE across 
CCD chips. After Source Extractor catalogs have 
been made, stars are picked out for their member¬ 
ship in the respective stellar loci of size-magnitude and 
magnitude-surface brightness diagrams. Initial identi¬ 
fication proceeds automatically using an algorithm that 
identifies the typical stellar size, then selects the mag¬ 
nitude range for which there is a significant density 
enhancement at that size. To ensure an accurate sam¬ 
pling of the PSE for circularization, the r band stellar 
catalogs are individually inspected and manually ad¬ 
justed as needed. The stars used for r band PSE fits 
are highlighted in the size-magnitude diagram of Eig- 
ure There and elsewhere in the paper, size is de- 
fineoas the sum of the second-order intensity moments 


























Mass Substructure in Abell 3128 


5 



Figure 1. Whisker plot showing spatial variation of the PSF across the Abell 3128 image. Each stick represents a star that was used to circularize the PSF, 
with length proportional to the magnitude of its measured ellipticity, and orientation equal to its position angle. Blank regions correspond to CCD chip edges, 
large cluster galaxies or saturated stars. Left: Stars in a Abell 3128 stack made without circularization correction. The mean segment size corresponds to stellar 
ellipticities ~ 0.05; points represent objects with e < 0.008. Right: Stack made with the multi-chip circularization correction described below. Mean ellipticity 
has been reduced to 0.005. 






Figure 2. Correlation functions between tangential ellipticity components for objects in the A3128 observation. Star-star auto-correlation functions are plotted 
at top left, and star-galaxy cross-correlation functions are contrasted with galaxy-galaxy auto-correlation functions in the other three panels. All ellipticity 
components are defined with respect to the image axes. 


hx + lyy wit h an additional factor to correct for non- 


Gaussianity ( Bernstein & Jarvis|2002 ) 


Master Catalogs. - The DLS survey found that the MS- 
CRED astrometric calibration of images is not good 
enough to stack them directly; small shifts between 
overlapping exposures would lead to spurio us stretch¬ 
ing of galaxy shapes ( [Wittman et al.|[2006| ). To pre¬ 
cisely define the astrometry of the final stack image, 


we match all the catalogs in equatorial coordinates to 
produce a master catalog. Every object that was ob¬ 
served in at least three exposures (within the tolerance 
of 1".8) has its mean right ascension, declination and 
magnitude recorded. Subsequently, the master catalog 
positions are used to define a coordinate system for the 
stack (a simple tangent plane projection with no optical 
distortion) and then transformed to pixel coordinates in 
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the final stack image. 

4. Pixel Coordinate Transformations. - For every CCD in 
each exposure, matches between the Source Extractor 
catalog and the master catalog are used to define a trans¬ 
formation from CCD pixel coordinates to final stack 
image. For the ugz^ data not subject to WL analy¬ 
sis, less stringent astrometry requirements allow for the 
DLS default of a third-order polynomial. In the case of 
r band exposures, the coordinate transformations must 
be defined by a fourth-order polynomial. Lower or¬ 
der polynomials underfit the variation at the edges of 
the CCD, resulting in a slight elongation of galaxies in 
the final stack image and ultimately leading to bands 
of spurious shear signal in the WL convergence maps. 
However, many CCDs from the edge of the exposure 
(~ 20%) have too few matches with the master catalog 
to support a 4th order fit and must be excluded from the 
final stack image. 

5. PSF Circularization. - Once the pixel transformation 
polynomials are determined, we generate PSF circular¬ 
ization profiles for every CCD image in the dataset. By 
default, the adaptive moment combinations Cx - fy, 
hx + lyy and Cy are fit automatically with 4th order 
polynomials. The exceptions are PSF profiles for the 
r-band imaging, which are prepared in adva nce us ing 
the multi-stellar catalog procedure laid out in |4.1.2 



r Magnitude 


Figure 3. Size-magnitude diagram of the object catalog generated from the 
final stacked r image. At this stage, the catalog is filtered for objects with 
Source Extractor or ELLIPTO error fiags, but no other cuts are applied. The 
stellar locus is the stripe of objects with size ~ 5 and 15.5 < nir < 21; red 
stars mark objects used in our multi-catalog PSF circularization procedure. 
Objects in the “second stellar locus” with sizes around 2 pixels and < 23 
are pixel noise variations. Spurious detectections along the bleed trails of 
saturated stars cause the uptick at m ~ 16. 

6. Photometric Calibration. - Before CCD exposures are 
stacked, their photometry is calibrated to ensure consis¬ 
tent object magnitudes everywhere in the final stack im¬ 
age. This is made more difficult since, due to the curva¬ 
ture of the focal plane, pixels at the edges of the DECam 
CCD subtend more sky area than pixels at the center of 
the array. To correct for focal plane distortion, catalog 
magnitudes are multiplied by the Jacobian of the pixel 
coordinate transformations computed in step 4 and then 
gathered into a master photometric catalog. For each 


CCD exposure, we then derive a relative photometric 
offset by matching its catalog to the master photomet¬ 
ric catalog and computing the 3 a clipped mean of the 
magnitude differences of the matching objects. 

Finally, to produce the stacked image, we implement the Deep 
Lens Survey’s DLSCOMBINE algorithm. For each pixel 
in the output image, DLSCOMBINE loops over contribut¬ 
ing pixels in the input images and applies the relevant bad 
pixel masks, PSF circularization kernels, coordinate transfor¬ 
mations and photometric offsets. A 3 cr clipping is applied 
before the mean pixel value is returned. A 3-color composite 
image made from the z, r and g stacked images is shown in 
Figure]^ 

4.2. Source Selection 

We produce a source catalog from the final sta cked r image 
using Source Extractor ( [Bertin & Arnouts|1996| ). The Source 
Extractor detection significance and deblending thresholds are 
deliberately set to low values so that the faint background 
galaxies on which we perform WL analysis will be recog¬ 
nized. To clean out the accompanying multitude of junk 
Source Extractor detections, we perform a number of cuts 
to the initial object catalog. All detections with high Source 
Extractor and ELLIPTO error flag values may easily be fil¬ 
tered out, but ridding the catalog of pixel noise “detections” 
presents a special problem. Because their few counts are con¬ 
tained within a small isophotal area, SExtractor frequently as¬ 
signs them reasonable magnitudes, and measurement of their 
adaptive moments produces no error fiags. To filter out such 
pixel noise from the object catalog, we removed detections 
with ELLIPTO-determined fluxes of less than 100 counts. 
The cuts on error fiags decreased the initial object catalog of 
1.4 million objects by about 22%, and the flux cut decreased 
it by a further 43%. 

The object catalog is then subject to the the size and magni¬ 
tude cuts typical of weak lensing studies, which must remove 
stars, low-redshift galaxies and noisy sources while maintain¬ 
ing a large sample. The criteria for inclusion in the final sam¬ 
ple were isophotal r magnitudes between 17.2 and 24.6, and 
object s ize be tween 6.0 and 200 pixeP, where size was de¬ 
fined in |4.1.2[ The atypically generous upper limit of the size 
cut reflects the fact that A3128 is at very low redshift, and its 
“background” contains many large galaxies. Requiring that 
objects be detected in all four bandpasses de facto constitutes 
an additional catalog cut, eliminating 30,000 objects from the 
final catalog. At this point, the filtering has reduced the cata¬ 
log down to 200,000 sources in total (25 sources arcmin"^). 

Any galaxies in the foreground of A3128 will not be 
sheared by the cluster, and so their presence in the final source 
catalog dilutes the convergence measured on the lens. To 
filter out low redshift contaminants, we derived photometric 
redshifts using BPZ ( |Bemtez|[2000| ) with the standard HDF 
prior. The ugrz magnitudes of cataloged objects were submit¬ 
ted to the program, although since the u-band observations are 
rather shallower than the grz observations, the photometric 
redshifts are essentially three-point fits. We used the default 
CWWSB template set (E, Sbc, Scd, Irr, SB3, and SB2) with 
no modifications, and allowed three levels of interpolation be¬ 
tween neighboring templates. The range to be considered was 
restricted to 0.03 < zbpz < 3.0. To evaluate BPZ results on 
low-redshift galaxies, we identified 10 galaxies in the obser¬ 
vations with spectroscopic redshifts z ~ 0.06 and used BPZ 
to determine their photometric redshifts. This test calls at- 
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Figure 4. Sequence of zrg composite images showing progressively higher magnifications of the Abell 3128 field. Top: The full DECam 1.5°x 1.5°field of view 
of the cluster and surrounding region. Left: Close-up view of the central 32’ x 30’ of Abell 3128. R ight: 3’ x 3’ image showing the background cluster ACT-CL 
J0330-5227, which hosts SUMSS J033057-52281, a radio source at z = 0.44 (jWerner et al.U2007J and references therein). The image is centered on a strong 
gravitational leasing arc associated with this distant cluster. 



Background Galaxy Redshift (z) 

Figure 5. Redshift distribution of the final catalog, filtered for low-redshift 
galaxies. The majority of our sources are at z ~ 0.30, with secondary peaks 
at z ~ 0.46 and z ~ 0.60. 

tention to the uncertainty in BPZ results, as all 10 galaxies 
were assigned redshifts between 0.07 and 0.13. This range is 
comparable to the per-galaxy rms error found by other studies 


(e.g. [Sehgal et al.||2008b . While most galaxies that are truly 
in the foreground of A3128 are large enough to be eliminated 
by the size cuts described above, we nonetheless took the low- 
redshift uncertainty into account when developing our redshift 
selection criterion of Zbpz > 0-19 at greater than 95% prob¬ 
ability. About 169,000 objects (21 sources arcmin"^) remain 
after this latest cut is made. 

One last step remains before we can make convergence 
maps for Abell 3128. Even after circularization, atmospheric 
seeing and the circularization procedure itself still smear 
out the measured adaptive moments, making galaxies appear 
rounder than they really are and diluting the WL shear signal. 
As a correction, each galaxy’s ellipticity is divided by a factor 
R wh ich relates the size of the g alaxy to the mean stellar PSF 
size ( [Bernstein & Jarvi^|2002| ). Galaxies must have a mea¬ 
sured R > 0.2 for the correction to be successful, otherwise 
they are cut from the catalog. Only about 5% of galaxies re¬ 
maining in the catalog at this stage fail to meet that criterion, 
which leaves a final source density of 20 sources arcmin"^. 

4.3. Convergence Mapping and Quantified Detection of 
Substructure 

The catalog finally contains the true shapes of confirmed 
galaxies behind A3128, and may be used to reconstruct the 
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Figure 6. An aperture mass peak in five maps of consecutively larger filter 
radii. Top: Detection significance of aperture mass in units of cr. Bottom: 
Signal-to-noise of the aperture mass. In this example adapted from our A3128 
observation, the signal peaks in both significance and signal-to-noise at a 
Schirmer hlter radius of 6500 pixels. 


cluster’s projected mass. To extract the aperture mass sig¬ 
nal of A3128 from the tangential ellipticities of background 
sources, this study relies on soft ware develope d by co-author 
Huwe as part of his thesis work ( Huwe||2013[ Huwe & Del- 
1 Antonio 2014, in preparation). Our particular implementa¬ 
tion of this software is presented here. 

Our software first bins the A3128 image into blocks 200 
pixels on a side to reduce computation time. Using the adap¬ 
tive moments of galaxies in the catalog, the software then re¬ 
turns the filtered aperture mass statistic (Equations and 
within each 200-pixel block. To quantify our mass recon¬ 
struction of the A3128 observation, we construct a WL signal- 
to-noise map as follows. Random noise maps are generated 
by recalculating the aperture mass statistic (with the same 
Schirmer filter and pixel block sizes) on a catalog of shuf- 
fied galaxy positions and moments. The random maps will 
initially be dominated by non-Gaussian shape measurement 
error, so the randomization process is repeated 100 times. As¬ 
suming that the errors in the aperture mass reconstruction will 
then be Gaussian, the variance of each image block in the ran¬ 
dom maps represents the 1 cr noise level of the Map statistic. 
Dividing the Map signal map by the variance of the random 
maps, pixel block by pixel block, yields an estimate for the 
signal-to-noise. 

Because WL distortion is tangential to the direction to the 
center of mass, an image should have no systematic B-mode 
(curl-like) distortion in the shapes of background galaxy. The 
lensing signal should thus vanish when Ctan, the E-mode (curl- 
free) component of shear, is replaced with the B-mode com¬ 
ponent Cc defined in Equation Any significant WL peaks 
obtained when ^tan is replaced with ec in Equation would 
not come from the cluster, but instead indicate some system¬ 
atic error in the analysis. Since most sy stematics are expe cted 
to add equal power to E- and B-modes ( Jarvis et al.|2003| ), we 
generate B-mode signal-to-noise maps to control for bias in 
our analysis. 

Both E- and B-mode signal-to-noise maps treat the errors 
in the convergence field as Gaussian, but this assumption may 
not be warranted. To further quantify confidence in the results 
of our mass reconstruction, we calculate the detection signif¬ 
icance of features in the signal-to-noise map with the follow¬ 
ing algorithm. Using the same filtered aperture mass statistic 
as before (Equation |^, the software creates a signal file and 
then iterates through some large number of random noise re¬ 
constructions which are stored in memory. At every 200-pixel 
block of the observation, the software tallies how many noise 


reconstructions had a greater magnitude of WL signal than 
the signal file. This number should be close to 0 for blocks 
near the cluster center, but in massless regions of the observa¬ 
tion will be roughly 50% of the total number of random iter¬ 
ations. When inserted into an inverse cumulative distribution 
function, this number is converted into a Gaussian-type con¬ 
fidence (T which quantifies the significance of the shear signal 
in that pixel block. The maximum attainable cr will depend on 
the number of noise iterations; our software generates 100000 
random maps which corresponds to a maximum confidence of 
4.42 cr. 

In addition to a magnitude given by its cr value, the soft¬ 
ware assigns to each significance map pixel the sign of the 
corresponding pixel block in the original signal map. Hence, 
regions in the significance maps with negative values corre¬ 
spond to statistically significant underdensities in the mass 
distribution, while positive cr means an area of mass enhance¬ 
ment compared to the mean. 

To search for mass concentrations, the significance maps 
are thresholded above -r3.8(T, and potential substructure peaks 
are identified by inspection. Eor each group of contiguous im¬ 
age blocks with significance greater than -r3.8 cr, we follow 
the feature through a range of Schirmer filter scales. Rela¬ 
tively small Schirmer filter radii do not encompass all the 
shear signal from the feature, whose significance will sub¬ 
sequently be suppressed. As the filter radius increases, the 
significance of the detection increases before peaking at some 
Schirmer filter size which is then the characteristic scale of 
that substructure. Eurther filter expansions eventually lead to 
the merging of the substructure signal into the overall cluster 
signal. An example of this increase and decrease in signifi¬ 
cance is shown in Eigurej^ 


4.4. NFW Shear Profile Fitting 

Recalling that aperture mass maps return only the relative 
mass enhancements in an observation, we have written soft¬ 
ware which fits A3128 and its substructures with axisymmet- 
ric NEW weak lensing shear profiles to constrain their physi¬ 
cal masses. Different algorithms are employed to fit the obser¬ 
vation with single and multiple NEW halos, but in both cases 
the software first obtains the scaling factor (Equation 
for each galaxy in the catalog. Requisite angular diameter 
distances for are computed from the galaxies’ BPZ red- 
shifts. The software then varies the M 200 of an NEW halo and 
computes the corresponding r 2 oo under a Planck XVI cosmol¬ 
ogy. The halo’s concentration c i s obtained by inserting its 
M 200 into the empirical relation of [Bhattacharya et 31] ( |2013 ) 
for their full cluster sample. The r 2 oo and c parameters be- 
come the and dc which characterize an NEW halo’s mass 
distribution and hence its shear profile. 

To fit a single NEW mass to the primary WL peak of 
A 3128, as M 20 Q is varied the s oftware follows the prescription 
of [Wright & Brainerd] ( |2000| to compute the halo’s reduced 
shear at the location of every background galaxy. We find the 
best-fit Mioo by u sing the parabolic extrapolation method of 


Press et al.| ( |2007| ) to minimizeresiduals between the NEW 
halo’s shear profile and the ^tan measured on the image. 

Simultaneous fitting of NEW masses to multiple substruc¬ 
tures requires a different approach: the (tensor) reduced shear 
from multi ple NEW peaks does not a dd linearly, so the pre¬ 


scription oflWright & Brainerd ( 2000|) is not directly applica¬ 
ble. Instead, we use the fact that background galaxies expe¬ 
rience an individual displacement fi from each NEW halo. In 
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the weak tensing approximation, the displacements from mul¬ 
tiple NFW haloes add linearly: Ptot = ZA- The total shear 
and convergence at every point in the image can then be built 
from derivatives of the Jacob ian (dPtotIdO) using the formu¬ 
lae of |Golse & Kneib| ( |2002| ). Presupposing the locations of 
their centers have been established, we vary the M 200 of NFW 
haloes centered on each substructure, obtain the correspond¬ 
ing reduced shear at the location of every background galaxy, 
and minimize the profiles’ residuals against the galaxies’ 
tangential ellipticities. 

We emphasize that neither the single-peak nor multiple- 
peak mass estimates in this work result from fitting 1-D 
NFW shear profiles to azimuthally averaged galaxy elliptic¬ 
ities, though such an approach is common in the literature. 
Rather, our NFW fitting method uses the full positional in¬ 
formation of every galaxy in the catalog. Tangential shear 
profiles shown below (Figure are for illustrative purposes 
only. We also note that in all mass estimates, NFW shear pro¬ 
files are centered on the aperture mass peak’s highest signal- 
to-noise pixel in c onve rgence maps. However, due to our bin¬ 
ning scheme (cf. ^4.3| ), each convergence map pixel actually 
spans 200 pixels (53") on the observation. The resultant am¬ 
biguity in the identified center of a WL peak could bias mass 
estimates through a mis-computation of the g alaxy ellipticity; 
we investigate this potential centroid bias in |5.4[ 


5. RESULTS 

5.1. Detection of Primary Cluster Aperture Mass 

Applying the procedure laid out in |4.3| to the observation, 
we report the detection of the Abell 3128 weak lensing sig¬ 
nal at high significance. The top left of Figure |7ja) shows 
the primary aperture mass peak, which saturates our signifi¬ 
cance maps with cr > 4.42 at all Schirmer filters larger than 
6000 pixels. To identify the aperture size which best charac¬ 
terizes the cluster, we constructed signal-to-noise maps using 
Schirmer filters up through Rg = 14000 pixels. Within this 
range, the A3128 aperture mass achieves its peak S/N of 8.4 
at two distinct locations with an aperture size of Rg = 10000 
pixels. Since the Schirmer filter weight peaks sharply at 
~ O.lORg, the primary A3128 signal spans roughly 4.4’ on 
the image. 

From the location of the highest signal pixels in the S /N 
map of Figure [TJb), we might assign the primary weak lens¬ 
ing peak of Abell 3128 to coordinates a = 3^30^50^5, S = 
-52°31'15". However, the presence of a massive high- 
redshift cluster (visible in Figure]^ only 6’ from the A3128 
X-ray center confounds the loca tion of its center of mass. In¬ 
stead, we defer this question to ^5.2| where the WL signal of 
the background cluster is probed. 

In addition to the primary cluster peak, several other high-cr 
aperture masses appear in Figure the 5 /A ~ 7 clumps at 
the bottom left of Figure |7jb) also saturate significance maps 
with cr = 4.42. These clumps are better ch aract erized at small 
Schirmer filter sizes, as described in in ^5.3| Immediately 
below A3128, the WL maps show a highly significant void in 
the local dark matter distribution. 

A signal-to-noise map made with the B-mode statistic Cc 
(Equatio n controls for any systematic error in the galaxy 
catalog ( |0| ). Accordingly, the B-mode signal-to-noise map 
for the Rg = 10000 pixel aperture is shown in Figure |^c). 
This map shows a marked anti-correlation between B-mode 
and E-mode signal in the neighborhood of the cluster, which 
becomes even more noticeable at small aperture sizes like Fig- 


ure[^ We investigate by making E- and B-mode maps from 
the catalog of stars used to circularize the PSE. In the absence 
of a systematic error, WL reconstructions from a catalog of 
perfectly circular(ized) stars would yield signal-to-noise maps 
indistinguishable from random noise. Instead, both the E- and 
B-mode S/N maps of Eigurej^ display a negative signal near 
the cluster center. Given that systematic errors should add 
equal power to E- and B-modes, it is expected to see the neg¬ 
ative signal in both the B-mode S/N map of Eigurej^c) and 
the E- and B-mode S/N maps of Eigurej^ 

Table 1 

Centroids from Previous Studies 


ID 

a 

(J2000.0) 

6 

(J2000.0) 

Reference 

A3128 Optical Center 
A3128 X-ray Center 
X-ray “SW Peak” 
ACT-CL J0330-5227 

3^30'”43" 

3 / 129 ^ 40 ^ 

3/129^55. 

3^30'”57" 

-52°31'50" 

-52°28'50" 

-52°34'50" 

-52°28'14" 

Rose et al. (2002) 
Werner et al. (2007) 
Werner et al. (2007) 
Menanteau et al. (2010) 


Equation shows that negative WL signal is caused by a 
net tendency of objects in a region to be aligned radially out¬ 
wards. Eurthermore, the DECam PSE is known to have a 
strong radial component (cf. the “pincushion” of Eigure[^. 
Hence, this negative signal is attributable to an undercorrec¬ 
tion of the PSE “pincushion” near the bright galaxies of the 
cluster, likely due to gaps in stellar coverage. This interpre¬ 
tation is supported by the weakening of the effect at large 
Schirmer filter radii, where the inclusion of more galaxies at 
random orientations washes out any localized residual corre¬ 
lations in the PSE. 

We constrain the magnitude of this systematic effect as fol¬ 
lows. We amalgamate the star and galaxy convergence maps 
and tally up WL signal enclosed in a circle centered on the 
cluster. This value is then compared to the equivalent signal in 
the convergence map made from galaxies alone. We find that 
the combined stars/galaxies convergence map has 5% more 
power than the galaxies-only signal map, likely because the 
systematic undercorrection of the PSE subtracted away some 
of the cluster’s original convergence signal. Since WL shear 
and mass scale linearly, we expect that the NEW fits presented 
in Table [^underestimate the true mass by ~ 5%. 


5.2. High-Redshift Background Cluster 

Visible in the last panel of Eigurej^is an interloping high- 
redshift cluster of galaxies, complete with the blue arc that is 
the hallmark of strong gravitational lensing. Multiple studies 
have confirmed that this background “cluster behind a clus¬ 
ter,” identified in the literature as ACT-CL J0330-5227, lies at 
^ = 0.44 and is the source of the north eastern lobe of X-ray 
emission in A3128 ( Werner et al.|2Q07| . All cataloged galax¬ 
ies with z > 0.44 will bear lensing signal from both clusters, 
and so the higher-redshift cluster must be precisely character¬ 
ized to prevent confusion with substructure in A3128. 

To disambiguate the signals of A3128 and ACT-CL J0330- 
5227, we prepared two subsamples using BPZ redshifts 
and probabilities: an intermediate-redshift sample containing 
galaxies at 0.19 < z < 0.4 with greater than 75% probability, 
and a high-redshift sample containing galaxies at z > 0.443 
with greater than 75% probability. The intermediate-redshift 
catalog contains 75,774 galaxies at a mean redshift of z = 
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Figure 7. Reconstructed weak tensing convergence maps for Abell 3128, all with 52’7pixel resolution on the observation. Left: Significance map made using a 
Schirmer aperture of radius 10000 pixels (= 44’). The primary peak is at the top left, and saturates our significance at 4.42cr. Center: Signal-to-noise map and 
44’ Schirmer aperture, which yielded the maximal cluster SIN of 8.4 at two separate locations. Right: B mode S/N map for 10000 pixel aperture radius (primary 
peak filter size). 
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Figure 8. Reconstructed aperture mass maps with 52” resolution made from the catalog of stars used to circularize the DECam PSF . Left: Signal-to-noise map 
made with a Schirmer aperture of radius 10000 pixels (= 44’), the filter size of the A3128 peak. Right: B mode S/N map for 10000 pixel aperture size. Note that 
both maps have significant areas of negative WL signal, indicating a net tendency of background galaxies to be alined radially outwards.) 

authentically identify the A3128 barycenter. From the loca¬ 
tion of the highest a pixel in large-aperture significance maps, 
we report the primary weak lensing peak (and presumably 
barycenter) of the cluster at a = 3^30^16^5, d = -52°33'57". 
This WL center is offset by 4.7 arcminutes from the optical 
center of the galaxy distribution, but by 6.24 arcminutes from 
the published X-ray center. Instead, the WL potential center 
coincides more closely with the “southwest peak” of X-ray 
emission. Having established a location for ACT-CL J0330- 
5227 and the A3128 barycenter in our WL reconstructions, 
we may proceed in identification of substructures in A3128 
with more assurance. 

5.3. High Significance Substructure 

Our significance maximization procedure yields two sub¬ 
structures within A3128 proper, which are visible in the left- 
hand panel of FigureThe peak to the cluster’s southwest 
(A2) saturates our significance maps at cr = 4.42 at apertures 
with Schirmer filters larger than Rg = 3500 pixels. The north¬ 
east peak (Al) does not achieve its maximum significance un¬ 
til Rs = 4500 pixels, at which point it, too, saturates our sig¬ 
nificance maps. The two substructures merge into one large 
4.42 (T structure in all significance maps made with apertures 
6000 pixels and larger. However, the two substructures can 
be still distinguished in S IN maps up through an aperture size 
of 9000 pixels. Both aperture masses achieve their maximal 


0.31, and the high-redshift catalog contains 52,847 galaxies at 
a mean redshift of z = 0.68. The intermediate-redshift galax¬ 
ies experience WL distortion only from A3128, since they lie 
behind A3128 but in front of ACT-CL 10330-5227; the high- 
redshift galaxies are behind both clusters and experience dis¬ 
tortion from both. 

In Figure significance maps made with the two redshift 
subsamples are compared to the the full galaxy sample map. 
The X-ray position of cluster ACT-CL 10330-5227 and the 
“southwest peak” of A3128 X-ray emission are marked with 
black stars. Using the high-redshift galaxy subsample (left 
panel), we detect at high confidence the WL signal from ACT- 
CL 10330-5227 and both the Al and A2 substructures of 
A3128. In the reconstruction with the intermediate-redshift 
galaxy subsample (center panel), the WL signal of A3128 is 
still distinct but the high-redshift cluster has dropped out of 
view. A reconstruction made with the full galaxy sample is 
shown in the right panel of Figure!^ the inclusion of a large 
number of galaxies at z < 0.44 dilutes the signal of ACT- 
CL 10330-5227 and makes it harder to discern. The fact that 
the background cluster does appear or disappear with such se¬ 
lections bolsters our confidence in the galaxies’ photometric 
redshifts, and by extension the angular diameter distances re¬ 
quired to fit NFW profiles to aperture masses. 

Based on the absence of high-z cluster signal in Figurej^b), 
the intermediate-redshift galaxy subsample may be us^ to 
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Figure 9. Significance maps marked with the published positions of the Abell 3128 X-ray “southwest peak” (bottom right star) and the z = 0.44 cluster ACT-CL 
J0330-5227 (top left star). Pixels in all maps span 52” on t he ob servation. The two extended high-cr regions apparent in all three panels are the principal 
substructures of A3128; they are discussed as A1 and A2 in g5.3| Left: Close-up of map made with background galaxy redshift restricted to z > 0.44 and a 
Schirmer aperture of 4000 pixels. The high-redshift cluster is plainly visible. Center: Map made with galaxies at redshifts between 0.16 < z < 0.4 and a 6000 
pixel aperture. While A3128 is still identifiable, albeit at slightly reduced significance, the high-redshift cluster has dropped out of view. Right: Map made with 
the full background galaxy sample and a 4000 pixel aperture. 


signal-to-noise at Rg = 8000 pixels: the southwest peak with 
SIN = 6.6, and the northeast peak with S/N = 5.6. Recall¬ 
ing that the Schirmer filter’s weight peaks at 0.1 ORy, the cor¬ 
responding angular scale of the substructures is 3.5 arcmin- 
utes. The northeast peak (Al) is spatially coincident with the 
brightest cluster galaxy, and both peaks are near the optical 
center of the cluster. Note that both of these features are spa¬ 
tially coincident with small-scale features in the “southwest 
peak” of X-ray emission (cf. Tables and 0, and are dis¬ 
tinct from the northeastern lobe of X-ray emission which is 
the high-redshift cluster ACT-CL J0330-5227. 

Within the same range of Schirmer apertures, we detect at 
high significance two substructures adjacent to the primary 
A3128 peak, circled in white in Figure [TT] and with coordi¬ 
nates listed in Table At Ry = 5500 pixels (corresponding 
to an object size of 2.4’), the topmost aperture mass (Gl) 
achieves its peak significance of 4A2cr while its neighbor¬ 
ing feature (G2) achieves AAlcr with a 5000 pixel aperture. 
At their respective characteristic aperture sizes, the features 
achieve S jN of 5.1 and 4.3. By Ry = 6000, both features 
have merged into the larger signal of A3128. From their spa¬ 
tial coi ncidence with knot s of galaxies at the same redshift as 
A3128 ( [Rose et al.|2002[ ), these two peaks are likely associ¬ 
ated with the cores of groups that have fallen into the clus¬ 
ter. 


Table 2 

Substructures Identified in Convergence Maps 


ID 

a 

(J2000.0) 

(J2000.0) 

Al/“Northeastern Peak” 

3'‘30"’02".6 

-52°38'33" 

A2/“Southwestern Peak” 

3'‘30"’28".4 

-52°33'57" 

B1 

3'‘35"'32".9 

-53°38T5" 

B2 

3'‘34"'50".0 

-53°35'40" 

B3 

3'‘32"’42".3 

-53°29'38" 

Gl 

3'‘3r'15".9 

-52°38'26" 

G2 

3*29»>55^o 

-52°34'50" 


The three high-cr features towards the bottom left of Fig¬ 
ure [TOTa) match up to regions of mass enhancement in the 
sign^to-noise image of Figurep^c). The leftmost peak (B1) 


is just at the edge of the detector, but peak B2 is less than an 
arcminute from the rich group AGO S 366 (z ~ 0.06). We 
surmise that peak B1 is caused by a combination of detector 
edge effects and a genuine DM enhancement associated with 
AGO S 366. Peak B1 achieves its maximum cr = 4.42 with a 
Schirmer filter size of 5000 pixels, and peak B2 does the same 
at Rs = 4500 pixels. The corresponding angular scales of 
these two aperture masses are 2.2’ and 2.0’, respectively. Be¬ 
ginning at peak B3, an arc of WL signal with 5/^-2.5-3.0 
stretches across the bottom of Figure [^b) and ends at the lo- 
cation of A3125 {a = 3''27.4'”,(5 = -^°31'). While there is a 
corresponding bridge of galaxies in the DECam image, there 
are only moderate enhancements in the significance map. We 
find no significant WL features at the location of A3125 itself. 

5 . 4 . Mass Estimates 

Having constrained the locations of the A3128 primary 
aperture mass and its substructures with WL convergence 
maps, we proceed to parametric ally fit them with NEW 
masses following the procedure in ^4.4[ 

5 . 4 . 1 . Single NEW fit to Primary A3128 Aperture Mass 

Centering a single NEW peak at the A3128 barycenter 
yields a best-fit mass of M 200 = 10.0 + 2.3 x The cor¬ 

responding r 2 oo for this mass is 2.1 Mpc, which at the distance 
of A3128 spans about 7000 pixels on the observation. Uncer¬ 
tainty is estimated through a jackknife approach wherein we 
randomly resample 50% of the galaxy catalog and recompute 
the best-fit mass. The variance of 2000 realizations is taken as 
a measure of the NFW fitting procedure’s internal consistency 
and becomes the error bar on the best-fit mass. 

5 . 4 . 2 . Simultaneous NFW fits to Abell 3128 and ACT-CL 

J0330-5227 

NFW shear profiles were fit to the A3128 barycenter 
and ACT-CL J0330-5227 simultaneously, using both the full 
galaxy sample and the high-redshift subsample {zgai > 0.44) 
from |5.2[ Parametric masses resulting from those fits are 
displayed in the left column of Table [funder the Centroid 1 
heading. We also performed fits with the A3128 shear pro¬ 
file centered on the southwest substructure peak (“Centroid 
2,” right column of Table [^. Errors on these and all other 
NFW fits in Table are omained through the variance of a 
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Figure 10. High resolution convergence maps made with with 52’7pixel resolution on the A3128 observation. Top Left: Significance map made with a Schirmer 
filter of Rs = 4500 pixels (= 19'), in which the two substructures attain 4.42 cr detection significance. Several peaks to the bottom left (Bl, B2 and B3) also 
attain high significance. Top Right: Magnified view of significance map at left, highlighting central cluster. Bottom left: E-mode signal-to-noise map made with 
a Schirmer filter of Rs = 4500 pixels. The Al and A2 substructure peaks are detected at *S/A > 5.6. Bottom Right: B mode S IN map made with Rs = 4500 pixel 
aperture. 



- 3.4 - 2.5 - 1.7 - 0.79 0.081 0.95 1.8 2.7 3.6 


Figure 11. WL significance map made with a Schirmer aperture of Rs = 5500 pixels (= 24'). At this filter size size, the two substructures in A3128 have nearly 
merged, but other structures surrounding the cluster are now visible. Inset: Magnified view of the WL significance map, overlaid on a zrg composite image of 
A3128. The two high significance peaks circled in white are likely associated with galaxy groups recently accreted onto the cluster. 


50% random resampling of the catalog. All fits were subject 
to 2000 resamplings except the full galaxy sample/Centroid 1 
fits, which were resampled 4000 times. 

In fits made with the full galaxy sample, the masses of the 
two clusters sum to about 2.2 x lO^^M©, but the allocation 
of this mass between A3128 and ACT-CL J0330-5227 varies 
significantly (> 1 cr error bars) depending on the A3128 cen¬ 
ter chosen. When the A3128 shear profile is centered at the 
barycenter, the high-z cluster is assigned a lower mass than 
when the A3128 profile is centered at the southwest substruc¬ 


ture peak. A similar picture emerges in the 2-peak fits with 
the high-redshift galaxy sample {zgai > 0.44): the combined 
mass of A3128 and ACT-CL J0330-5227 remains the same 
(there, 2.7x lO^^M©) regardless of A3128 centroid chosen, but 
the distribution of this mass between the two clusters varies 
significantly. We note that compared to the fits with the full 
galaxy sample, the variance of the A3128 and ACT-CL J0330- 
5227 masses in Table are slightly higher - an expected by¬ 
product of the smaller sample size. 

Table also contains the results of a two-peak fit made 
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with the inter media te-redshift subsample of galaxies (0.16 < 
Zgai < 0-4) of ^5.2[ As expected when using galaxies exclu¬ 
sively in front ot the high-redshift cluster, no mass was allo¬ 
cated to ACT-CL J0330-5227. The mass assigned to A3128 
in this round of two-peak fits underestimates the other val¬ 


ues in Table at the > 1.5 cr level. Because of the low 
source density of the intermediate-redshift subsample (7-10 
arcmin"^), and moreover because those galaxies have an unfa¬ 
vorable Ds /Dls , the underestimate likely reflects a WL signal 
too weak to beat down the random shape noise of background 
galaxies. 


Table 3 

NFW Masses for High-Significance Aperture Masses 



A3128 Centroid 1 

A3128 Centroid 2 

Galaxy Sample 

A3128 Mass^* 
(lO'-’A/o) 

High-z Mass 
(10‘‘*Mo) 

A3128 Mass'’ 
(1O“'M0) 

High-z Mass 
(lO'-’A/o) 

Full Zgai sample 
Zgal ^ 0.44 

0.16 < Zgal < 0.4 

11.5 + 1.7 
13.8 + 2.4 
7.0+ 1.9 

10.3 + 3.2 

14 + 4.0 
<0.4^ 

7.3 + 1.3 

8 .0+ 2.1 

3.5 + 1.3 

15.3 + 2.4 
19.5 + 3.5 
<0.4 


^ Centered at cluster barycenter: a = 3^30"^ W.5,6 = -52°33'57" 

^ Centered at SW substructure peak: u = 3^30'”02^6, 6 = -52°38'33" 

Fits with this galaxy sample consistently assigned to the high-redshift cluster mini¬ 


mum mass in the allowed range, down to 1 x 10^^M© 


The;^^ landscape shown in Figure 12 a) indicates that the 
A3128 mass is well constrained by ^ full galaxy sample 
fits, but that they only place marginal constraints on the high- 
redshift cluster mass. Eliminating all galaxies in front of the 
high-redshift cluster leads to slightly better constraints on the 
mass of ACT-CL J0330-5227, as evidenced by the tighter 
ellipse about the best-fit masses in Figurep^b). Coupled with 
the fact that the high-z cluster mass experiences significantly 
greater variance in all fits. Figure [T^ suggest that A3128 dom¬ 
inates the WL signal embedded iiTthe background galaxy cat¬ 
alogs. 


5 . 4 . 3 . NFW Fits to Substructures 


We also attempted to fit masses to the two high- 
significance, small-aperture substructures identified within 
the A3128 peak (A1 and A2 in Table [^. The lowest fit 
assigned a mass of L02x lO^^M© to the northeastern substruc¬ 
ture and the minimum boundary value mass of 0.5 x IO^^Mq 
to the southwestern. These results were replicated in a simul¬ 
taneous 3 NFW shear profile fit which also included ACT-CL 
J0330-5227. Although the combined masses of the two sub¬ 
structures matches the best-fit A3128 mass in Table their 
allocation appears inconsistent with the larger size and higher 
S /N of the southwestern substructure and the fact that it satu¬ 
rates significance maps sooner than the northeastern substruc¬ 
ture. The shunting of the A3128 mass to the northeast aper¬ 
ture mass likely reflects its proximity to the cluster barycenter 
(400 pixels away) compared to the southwest aperture mass 
(over 1000 pixels away). Figure p^c) supports an interpre¬ 
tation where the total mass of the cluster is constrained to 
1.0 X 10^^ Mq by this set of fits, but that the two substructures 
cannot be resolved with individual NFW profiles. 

Finally, NFW masses were fit to the two infalling groups 
identified in Figure pT] (G1 and G2 in Table in a simul¬ 
taneous 4 NFW shear profile fit with A3128 and ACT-CL 
J0330-5227. To guarantee the distinctness of G1 and G2 from 
the shear profiles of the two clusters, A3128 and ACT-CL 
J0330-5227 masses were fixed to their respective full sam- 
ple/Centroid 1 values of 1.1 x IO^^Mq and 1.03 x lO^^M©. 
The fits returned masses of 2.7 + 3.2 x lO^^M© for Gl, and 
2.2 + 1.7 X IO^^Mq for G2, where the error bars come from 


the variance of 2000 jackknife resamplings, as before. 

5 . 4 . 4 . Global Tests of NFW Fitting Procedure 

In Figure[^ the best-fit A3128 single-peak tangential shear 
profile is compared with the azimuthally averaged galaxy el- 
lipticity signal. We again emphasize that mass estimates in 
this work do not result from fitting 1-D NFW profiles to 
binned galaxy ellipticities; the shear profiles shown below are 
for comparison only. The wide area of the A3128 observation 
and large number of background galaxies allows for a fine ra¬ 
dial binning and detailed inspection of galaxy ellipticity sig¬ 
nal. The negative value in the first radial bin is a manifestation 
of t he P SF undercorrection near the cluster center discussed 
in |5.1| At large R - Rc, the galaxy ellipticity signal should 
approach zero; however Figure p3] sho ws a noticeable down¬ 
ward trend in galaxy ellipticity. As in |5.1[ we attribute this 
negative galaxy signal to an undercorrection of the PSF “pin¬ 
cushion,” manifesting as a net tendency of galaxies at large 
distances from the cluster center to be aligned radially out¬ 
wards. To avoid underestimating the WL signal in A3128 
when computing the masses shown in Table we truncated 
the catalog to galaxies 58 arcminutes from the cluster center. 
Compared to NFW fits using an untruncated galaxy sample, 
this had the effect of raising the masses by 10%. 

Uncertainties on cluster masses in Table [^presume that our 
measurements obey Gaussian statistics. This assumption may 
be tested using the distribution of cluster masses returned by 
the jackknife procedure. The left panel of Figure p^ shows 
a 2-D histogram of cluster masses returned by 4000 random 
resamplings of the full galaxy catalog. The right panel of Fig¬ 
ure [14] shows the distribution of allowed A3128 NFW masses 
returned in the resampling, i.e., the left panel collapsed along 
the y-direction of ACT-CL J0330-5227 mass. Starting at the 
the median mass (~ l.Ox lO^^M©), we sum 34.1%, 47.7% and 
49.9% of the returned masses on either side of the distribution. 
The equivalent 68% confidence interval is (7.6,12.2)x 10 ^"^Mq 
and the 95% confidence interval is (5.9,15.6)x 10 ^"^Mq. When 
summing over 99% of the returned mass range, the bound on 
the low-mass end is 4.6x lO^'^M©. However, the upper bound¬ 
ary value of the sampled mass range is reached before we can 
find an equivalent upper bound. The 68% confidence interval 
is symmetric about the median (i.e., LO^^^ x lO^^M©), and is 
also roughly equivalent to the 1 cr variance of Table How¬ 
ever, the 95% and 99% confidence intervals are not equiva¬ 
lent to 2 and 3 cr, nor are they symmetric about the median. 
This skewness to high masses suggests that the distribution 
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Figure 12. landscapes for fits of NFW shear profiles to the WL peaks identified in the reconstructions. For ease of viewing, the values have been rescaled to 
each distribution’s respective (x^ ~Xmin^ ^ 1000. Left: Residuals from parametric mass fits to ACT-CL J0330-5227 and A3128 (centered at its barycenter) made 
with the full galaxy sample. Center: Residuals from parametric mass fits to ACT-CL J0330-5227 and A3128 (centered at its barycenter) with Zgai > 0.44. Right: 
Residuals from fits to the two A3128 substructures with the mass of ACT-CL J0330-5227 fixed to 1.0 x lO^^M©. 



Figure 13. Tangential shear profile for the 1-peak NFW fit to the A3128 
barycenter (solid red line), overplotted on the azimuthually averaged tangen¬ 
tial ellipticity signal of background galaxies (solid blue line). The dashed line 
is the B-mode signal of galaxy ellipticity. Error bars on the galaxy ellipticity 
signal are the value of reduced shear in a radial bin (~ <^tan>/ V2) divided by 
yfN galaxies in that bin. The best-fit mass of 1.1 x lO^^M© was used for the 
theoretical curve. 

of allowed masses departs from Gaussianity at the endpoints, 
though it is Gaussian near the best-fit value. 

The finite resolution of our significance maps leads to to an 
uncertainty in the coordinates of WL peak centroids. Any re¬ 
sulting mis-centering of NFW profiles might bias the reported 
masses. This potential systematic was probed through a slew 
of mass fits in which the identified centers of A3128 and ACT- 
CL J0330-5227 were individually shifted north, south, east, 
and west by 200 pixels (the size of the centroid uncertainty). 
The mass of ACT-CL J0330-522 varied an average of 7.3% 
from its Table value, while the mass of A3128 only varied 
by 2.7%. Since the variations are smaller than the error bars 
reported in Table the uncertainty of centroid coordinates is 
probably not a dominant source of error in our results. 

6 . DISCUSSION 

With increasing Rg, the smaller aperture mass peaks de¬ 
tected in the periphery of A3128 such as ACO S 366 (B1 


in Figure display the gradual in creas e and decrease of 
significance and S IN expected from However, the pri¬ 
mary cluster saturates significance maps with cr = 4.42 at all 
Schirmer radii considered. In addition, the two A3128 sub¬ 
structures merge into the primary cluster signal at Rg = 6000 
in significance maps, but remain distinct in S /N maps through 
much larger kernel sizes. These results suggest that the pri¬ 
mary cluster and its substructures are detected at much higher 
confidence than 4.42 cr. They also underscore a fundamental 
limitation of the significance maps, which are only as good 
as the number of random noise iterations performed. Un¬ 
fortunately, the gain in significance is a slow function of the 
number of random maps. For example, to achieve cr = 4.89 


requires 1 x 10^ iterations, and cr = 5.33 requires 1 x 10^ 
iterations or 100 times more computation time. Given unlim¬ 
ited computational resources, the two substructures could be 
teased apart from the primary cluster signal at Rg > 5000 
pixels and firmly assigned an equivalent cr of detection confi¬ 
dence. 

We note that the departure from Gaussianity in Figure 
does not invalidate our use of significance maps to evaluate 
WL signal. The “cr” of a significance map is a Gaussian- 
equivalent confidence, representing a signal pixel’s distance 
from the mean , and is rooted in an exact pixel probability dis¬ 
tribution (see §4.3| ). Should the distribution of random noise 
pixels have the same skew positive seen in Figure the cr 
of our significance maps - as well as our S /N values, which 
assume a Gaussian noise distribution - would in fact only 
be pseudo-Gaussian. However, the relative pixel-to-pixel en¬ 
hancements in significance maps would still be genuine: a 2 
cr peak is still 95.4% less likely than the mean value of WL 
signal, a 3 cr peak is still 99.7% less likely than the mean value 
of WL signal, etc. For this reason, the significance maps and 
the distribution of allowed NFW masses are complementary: 
while jackknife resampling makes no assumptions about the 
distribution of WL signal pixels, it must be remade for each 
aperture mass as it contains no spatial information. On the 
other hand, the significance maps (whether strictly Gaussian 
or not) do reveal the 2-D distribution of projected mass. 


While ACT-CL J0330-5227, A3128, and the A3128 sub¬ 


structures are all unambiguously detected in WL reconstruc¬ 
tions at very high significance, their parametric mass fits are 
subject to some degeneracy. Both in simultaneous NFW fits 
to the two clusters and in simultaneous fits to the A3128 sub- 
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\J.^ \J.\J \J.u ±.\J ±.Z. u.^ U.U U.O l.U l.Z 1.^ l.U 

A3128 Mass (IO'^A/q) A3128 Mass (lo'^A/g) 

Figure 14. Left: 2-D histogram of parametric masses for A3128 (with NFW shear profile centered at the barycenter) and ACT-CL J0330-5227 returned by 4000 
jackknife resamplings of the full (untruncated) galaxy catalog. Right: Distribution of most probable NFW masses for A3128 with an unconstrained ACT-CL 
J0330-5227 mass. This plot is equivalent to projecting the 2-D histogram at left along its y-axis of ACT-CL J0330-5227 masses. The spike at the end of the 
x-axis is caused by projecting all masses greater than 1.55 x 10^^ into the highest bin. 
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structures, the sum of the masses is more tightly constrained 
than their respective magnitudes. This is evident from all 
three plots of Figure 12 in which the best-fit ellipses tilt 


about a line of constant mass. Table |3] also attests to the de¬ 
generacy of mass fits: with both the full and high-redshift 
galaxy samples, the aggregate mass of ACT-CL J0330-5227 
and A3128 is conserved, but their respective magnitudes vary 
at the < 1 (T level depending on the chosen A3128 NFW pro¬ 
file centroid. The significant dependence of cluster mass on 
A3128 centroid signals that the exact assignment of mass in 
our two-peak fits is driven by sources near the center of the 
cluster, and that galaxies further away only bear the overall 
signal. As it happens, the source density in our observation is 
depressed near the cluster center - a consequence of studying 
low redshift clusters whose member galaxies subtend wide 
angles. This explains the degeneracy in our parametric mass 
fits, and particularly our inability to resolve the A3128 sub¬ 
structures with NFW shear profiles. 

Figure shows a downward trend in galaxy ellipticity at 
distances past 50-60 arcminutes from the cluster center, which 
we attributed above to an under-correction of the “pincush¬ 
ion” in the DECam PSF. A residual radial gradient in the PSF 
was also invoked explain the non-zero correlation at small 
scales visible in Figure and the negative galaxy signal in 
Figure To avoid underestimating the WL signal in the ob¬ 
servation, the NFW mass fits were supplied exclusively with 
galaxies within 58 arcminutes of the cluster center. In fu¬ 
ture work, we will attempt to improve our PSF circulariza¬ 
tion sc heme, perhaps adopting the method of |Miyatake et al. 


( 2013| ) to reject exposures with unsatisfactory PSF correction. 

Aside from highlighting the shortcomings of our PSF cor¬ 
rection scheme, Figure[^has another interesting feature. The 
large number of galaxies in our observation allows for a fine 
binning of the galaxy ellipticity signal which, in turn, allows 
us to appraise the success of an NFW profile in describing it. 
Ignoring the inner 5 arcminutes, it is clear that the galaxy el¬ 
lipticity distribution tends to zero faster than predicted by an 
NFW fit: galaxy ellipticities are already consistent with zero 
by 35 arcminutes, while the NFW fit does not settle to zero 
on the scale of the imag e. Cosmological simulation s such as 
have been published by |Diemer & Kravtsov| (|2014|) report a 
similar conclusion, that the density profiles of cluster-size ha¬ 
los fall off more quickly than predicted by either the NFW or 


Einasto mass models. In future work, we will attempt to fit 
the galaxy ellipticity distributions with the latest generation 
of mass distributions. 


7. CONCLUSIONS 

The significance maps and the distribution of most probable 
masses in Figure [T^ are completely independent ways of con¬ 
firming the WL signal from Abell 3128. Significance maps 
employ the aperture mass statistic of Equation to sum up 
background galaxies’ tangential ellipticity and return the con¬ 
vergence signal at each point in the observation. Eigure 
results from the parametric fitting of an NEW halo shear pro¬ 
file to the tangential ellipticities of background galaxies. With 
either method, the weak lensing signal from Abell 3128 is de¬ 
tected at high significance. In particular, the probability dis¬ 
tribution of A3128 mass obtained by randomly resampling of 
the full galaxy catalog skews to high masses, from which it 
follows that low values of A3128 mass are more strongly dis¬ 
favored than higher masses. Our confidence in the WL de¬ 
tection of A3128 is thus higher than might be indicated by 
significance maps alone. 

Given a sufficient density of background galaxies (and with 
the caveat that the number of degrees of freedom decreases 
with added NEW peaks), the tools we have developed allow 
for the fitting of NEW masses to an arbitrary number of sub¬ 
structures. Our average source density of 20 galaxies per ar- 
cminute was enough to constrain the mass of Abell 3128 to 
(1.1 + 0.2) X IO^^Mq through both the simultaneous NEW 
fits with the high-redshift cluster ACT-CL J0330-5227 and the 
substructure mass fits. However, the background galaxy den¬ 
sity drops to ~ 13 arcmin"^ near the cluster center because 
of the large apparent size of its member galaxies. Combined 
with the unfavorable distance ratio of background galaxies 
in our observations, these issues prevent the fitting of NEW 
masses to the two central substructures. Moving forward, 
we expect that all WL studies of low-redshift galaxy clus¬ 
ters will be similarly affected, and this places a lower limit 
on the source density required to resolve their substructures 
with NEW shear profiles. 

Since they can be mistaken for mass substructures, inter¬ 
loping high-redshift clusters like ACT-CL J0330-5227 pose a 
problem for systematic studies of low-redshift clusters. Even 
when they are not fully resolved in WL convergence maps. 
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Figure 15. Abell 3128 raw convergence map for a 10000 pixel Schimer aperture size, colorized and superimposed on zrg composite image. The bright clump 
on the top left is Abell 3128; the fainter clumps on the bottom left are from infalling groups at z = 0.06 - 0.08. The dispersed arc of signal along the bottom is 
from Abell 3125, a cluster disturbed by a recent passage near Abell 3128 (Werner et al 2007). 


high-redshift background clusters can have knock-on effects 
on NFW mass estimates because their presence may blur the 
location of the low er-redshift cluster’s barycenter. The tests 
performed in ^5.2| were a means of vetting the extraction of 
photometric redshifts from deep DECam imaging, as we ex¬ 
pected to see WL signal from ACT-CL J0330-5227 only in 
galaxies with z > 0.44. The success of the redshift tests sug¬ 
gests a natural means of authenticating potential mass sub¬ 
structures: a localized WL signal enhancement that appears 
only when distant galaxies are used - like ACT-CL J0330- 
5227 in Figure [^ - is almost certainly an interloping high- 
redshift cluster. 

Abell 3128 is one of the lowest redshift clusters to have 
been studied with weak gravitational leasing in such detail. 
The advent of wide-angle cameras such as DECam makes the 
systematic studies of low-redshift clusters possible. In partic¬ 
ular, the work on A3128 presented here is the pilot for our WL 
study of the mass distributions of a complete, volume-limited 
sample of massive galaxy clusters between 0.04 < z < 0.1 
using DECam. We hope that this work will in turn become 


the low-redshift anchor to a systematic observational mea¬ 
surement of the evolution of mass substructure in clusters of 
galaxies. The first of its kind, such a study would deepen our 
understanding of how the first galaxies assembled themselves 
into clusters. In addition, the presence of substructure may in¬ 
duce biases in the mass determination of clusters that must be 
precisely calibrated to extract precision cosmological infor¬ 
mation from the cluster mass function and its evolution. Up¬ 
coming surveys such as the Dark Energy Survey (DES) and 
LSST will measure averaged lensing signals from thousands 
of clusters to calibrate the cluster mass function. Therefore, 
empirically measuring the role of substructure in mass mea¬ 
surements is now of great importance. 

The authors thank Dara Norman for her efforts in per¬ 
forming all observations made during the DECam science 
verification phase, including those of Abell 3128 used in 
this study. The authors also thank the anonymous referee 
for suggesting many useful improvements to the analysis. 
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Champaign, Institut de Ciencies de I’Espai, Institut de Eisica 
d’Altes Energies, Lawrence Berkeley National Lab, Ludwig- 
Maximilians Universitat, University of Michigan, National 
Optical Astronomy Observatory, University of Nottingham, 
Ohio State University, University of Pennsylvania, University 
of Portsmouth, SLAC National Lab, Stanford University, Uni¬ 
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lence “Origin and Structure of the Universe” and the DES col¬ 
laborating institutions. 

REFERENCES 


Andernach, H., Plionis, M., Lopez-Cruz, O., Tago, E., & Basilakos, S. 2005, 
Nearby Large-Scale Structures and the Zone of Avoidance, 329, 289 
Bartelmann, M., & Schneider, R 2001, Phys. Rep., 340, 291 
Benitez, N. 2000, ApJ, 536, 571 
Bernstein, G. M., & Jarvis, M. 2002, AJ, 123, 583 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 

Bhattacharya, S., Habib, S., Heitmann, K., & Vikhlinin, A. 2013, ApJ, 766, 
32 

Clowe, D., De Lucia, G., & King, L. 2004, MNRAS, 350, 1038 
Coe, D., Benitez, N., Broadhurst, T., & Moustakas, L. A. 2010, ApJ, 723, 
1678 

DePoy, D. L., Abbott, T., Annis, J., et al. 2008, Proc. SPIE, 7014, 

Diemer, B., & Kravtsov, A. V. 2014, arXiv: 1407.4730 


Duffy, A. R., Schaye, J., Kay, S. T., & Dalla Vecchia, C. 2008, MNRAS, 

390, L64 

Fahlman, G., Kaiser, N., Squires, G., & Woods, D. 1994, ApJ, 437, 56 
Fischer, R, & Tyson, J. A. 1997, AJ, 114, 14 

Flaugher, B. L., Abbott, T. M. C., Angstadt, R., et al. 2012, Proc. SPIE, 

8446, 

Fleenor, M. C., Johnston-Hollitt, M., Rose, J. A., Christiansen, W. A., & 
Hunstead, R. W. 2008, Frontiers of Astrophysics: A Celebration of 
NRAO’s 50th Anniversary, 395, 367 
Gao, L., White, S. D. M., Jenkins, A., Stoehr, R, & Springel, V. 2004, 
MNRAS, 355, 819 

Gao, L., Navarro, J. R, Frenk, C. S., et al. 2012, MNRAS, 425, 2169 
Golse, G., & Kneib, J.-P. 2002, A&A, 390, 821 
] Hetterscheidt, M., Erben, T., Schneider, R, et al. 2005, A&A, 442, 43 
Hodge, J. A., Karim, A., Smail, I., et al. 2013, ApJ, 768, 91 
Huwe, P. M. 2012, American Astronomical Society Meeting Abstracts #220, 
220, #435.10 

Huwe, P. M. 2013, Dark Matter Substructure in High Redshift Clusters of 
Galaxies, Ph.D. Thesis, Brown University 
Jarvis, M., Bernstein, G. M., Fischer, R, et al. 2003, AJ, 125, 1014 
Kaiser, N., & Squires, G. 1993, ApJ, 404, 441 
Kaiser, N., Squires, G., & Broadhurst, T. 1995, ApJ, 449, 460 
Klypin, A. A., Trujillo-Gomez, S., & Primack, J. 2011, ApJ, 740, 102 
Lau, E. T., Nagai, D., & Nelson, K. 2013, ApJ, 777, 151 
Massey, R., Refregier, A., Bacon, D. J., Ellis, R., & Brown, M. L. 2005, 
MNRAS, 359, 127 

Miyatake, H., Nishizawa, A. J., Takada, M., et al. 2013, MNRAS, 429, 3627 
Nelson, K., Lau, E. T., Nagai, D., Rudd, D. H., & Yu, L. 2014, ApJ, 782, 107 
Okabe, N., & Umetsu, K. 2008, PASJ, 60, 345 

Okabe, N., Futamase, T., Kajisawa, M., & Kuroshima, R. 2014, ApJ, 784, 90 
Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2013, 
arXiv: 1303.5076 

Plazas, A. A., Bernstein, G. M., & Sheldon, E. S. 2014, Journal of 
Instrumentation, 9, 4001 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 2007, 
Numerical recipes in C (3rd ed.; New York: Cambridge University Press) 
Rose, J. A., Gaba, A. E., Christiansen, W. A., et al. 2002, AJ, 123, 1216 
Schirmer, M., Erben, T., Schneider, R, Wolf, C., & Meisenheimer, K. 2004, 
A&A, 420, 75 

Schneider, P. 1996, MNRAS, 283, 837 

Sehgal, N., Hughes, J. R, Wittman, D., et al. 2008, ApJ, 673, 163 
Taniguchi, Y, Scoville, N., Murayama, T., et al. 2007, ApJS, 172, 9 
Werner, N., Churazov, E., Finoguenov, A., et al. 2007, A&A, 474, 707 
Wittman, D. 2002, in Lecture Notes in Physics, ed. Courbin, F. and Minniti, 
D. (Berlin: Springer-Verlink), 608, 55 
Wittman, D. M., Tyson, J. A., Dell’Antonio, I. R, et al. 2002, Proc. SPIE, 
4836, 73 

Wittman, D., DelTAntonio, I. P, Hughes, J. R, et al. 2006, ApJ, 643, 128 
Wright, C. O., & Brainerd, T. G. 2000, ApJ, 534, 34 


